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The pressure dependence of the magnetic penetration depth A in polycrystaUine samples of 
YBa2Cu30a; with different oxygen concentrations x = 6.45, 6.6, 6.8, and 6.98 was studied by muon 
spin rotation (/iSR). The pressure dependence of the superfluid density pg oc l/\^ as a function of 
the superconducting transition temperature Tc is found to deviate from the usual Uemura line. The 
ratio {dTc/dP)/{dps/dP) is factor of ~ 2 smaller than that of the Uemura relation. In underdoped 
samples, the zero temperature superconducting gap Ao and the BCS ratio Ao/fcsTc both increase 
with increasing external hydrostatic pressure, implying an increase of the coupling strength with 
pressure. The relation between the pressure effect and the oxygen isotope effect on A is also dis- 
cussed. In order to analyze reliably the /.iSR spectra of samples with strong magnetic moments in 
a pressure cell, a special model was developed and applied. 



PACS numbers: 



I. INTRODUCTION 



The compound YBa2Cu30i: was the first high tem- 
perature superconductor— (HTS) with a superconducting 
transition temperature Tc above the boihng point of liq- 
uid nitrogen, and is one of the most studied HTSs<^ Its 
superconducting properties arc well characterized, even 
though some of them arc still being heavily discussed. 
Detailed muon spin rotation (/iSR) studies of the mag- 
netic penetration depth A and the superfluid density 
Ps oc 1/A^ were performed on poly- and single crystals of 
YBa2Cu30a; at ambient pressurci^"— However, the key 
question concerning the pairing mechanism responsible 
for high temperature superconductivity is still not re- 
solved, and is subject of intense debates. Although it is 
widely believed that magnetic fluctuations play a domi- 
nant role in the pairing mechanism,— oxygen isotope ef- 
fect (OIE) studies indicate that lattice degrees of freedom 
are essential for the occurence of superdonductivityJ^"— 
By means of isotope substitution one can probe the in- 
fluence of lattice degrees of freedom on superconductiv- 
ity without changing the lattice parametersi^ There are 
no other easily accessible methods which allow to solely 
modify the exchange integral J, in order to investigate 
its influence on the superconducting statci^ However, 
the application of hydrostatic pressure changes the in- 
teratomic distances in the lattice which in turn modifies 
both the lattice dynamics^- and the exchange coupling J 
between the Cu spins in cuprates i^^i^^ Therefore, a de- 
tailed study of the pressure effect (PE) on the supercon- 
ducting properties, e.g., the superfluid density Ps oc 1/A^, 
the gap magnitude Aq, and the BCS ratio Ao/ksTc, may 
provide important information for testing microscopic 
theories of the high-temperature superconductivity! ^^'^^ 

Up to now, the PE on the superconducting transi- 
tion temperature Tc was studied by resistivity and Hall 
effect experimentsi^"— Several phenomcnological^Si^iS 



and microscopic models were proposed based on a 
Hubbard^i^ or a general BCS approach in order to ex- 
plain the PE on Tci^ The role of nonadiabatic effects 
is discussed in Ref. [S^- These models suggest two ba- 
sic sources for the PE on Tc'. (i) A charge transfer from 
the charge reservoir to the superconducting Cu02 plane, 
which was conflrmed by Hall effect experimentsj^i^ and 
(ii) an increase of Tc due to a pressure dependent pairing 
interaction. 

The magnetic penetration depth A is a fundamen- 
tal parameter of a superconductor. It is a measure of 
the superfluid density according to the relation 1/A^ oc 
ris/m*, where n^ is the superconducting carrier dan- 
sity and m* is the corresponding effective massi^ From 
the temperature or flcld dependence of A one can de- 
termine the symmetry of the superconducting gap. its 
magnitude and the BCS ratio. The pressure dependence 
of A was previously studied in fine powdered grains of 
YBasCuaO:^ and YBa2Cu40a^-ii by means of mag- 
netization experiments. The ^SR technique is powerful 
and direct method to determine A in the bulk of a type- 
II superconductor 42il^ However, due to several technical 
difficulties only a small amount of /xSR studies of the 
penetration depth under pressure were performed so far. 
The main technical problems are: (i) The low fraction of 
muons stopping in the sample inside the pressure cell and 
(ii) the strong diamagnetism of a superconductor which 
substantially influences the ^SR response of the pressure 
cell. 

Here, we report on pressure dependent magnetic pen- 
etration depth studies in polycrystaUine samples of 
YBa2Cu302; {x = 6.45, 6.6, 6.8, and 6.98) by means of 
/iSR. We found that the pressure-dependent superfluid 
density ps oc 1/A^ vs Tc does not follow the Uemura 
relation.i^ The ratio Up = {dTc/dP)/{dps/dP) is a fac- 
tor ~ 2 smaller than that of the Uemura relation, but is 
quite close to that found in oxygen isotope effect (OIE) 



studiesr^iii suggesting a strong influence of pressure on 
the lattice degrees of freedom. Interestingly, a small pres- 
sure dependence of the superluid density was also found 
in the overdoped sample {x = 6.98). The superconduct- 
ing gap Aq and the BCS ratio Ao/fesTc both increase 
upon increasing the hydrostatic pressure in the under- 
doped samples, hence implying an increase of the cou- 
pling strength with pressure. Finally, a method of data 
analysis for tranverse-field /iSR measurements of mag- 
netic/diamagnetic samples loaded in a pressure cell is 
presented and applied here. This method leads to a sub- 
stantial reduction of systematic errors in the data analy- 
sis. 

The paper is organized as follows: In Sec. II we give 
some experimental details. In Sec. Ill we describe the 
method of /iSR data analysis and present the experimen- 
tal results, followed by a discussion in Sec. IV. The con- 
clusions are given in Sec. V. In the Appendix we describe 
the method used in this work in order to analyze ^SR 
spectra obtained for a magnetic/superconducting sample 
loaded in a pressure cell. 



II. EXPERIMENTAL DETAILS 

High quality polycrystallinc YBa2Cu30^ samples with 
X = 6.98, 6.8, 6.6, and 6.45 were prepared from the 
starting oxides and carbonate Y2O3, CuO and BaCOs 
as described elsewhere4i Transverse field (TF) iiSR ex- 
periments were performed at the ^El and 7rM3 beam 
lines of the Paul Schcrrer Institute (Villigcn, Switzer- 
land). The samples were cooled in TF down to 3K, and 
/iSR spectra were taken with increasing temperature in 



applied fields B^ 



0.1 and 0.5 T. Typical statistics 



for a ^SR spectrum were 5 — 6 x 10 positron events in 
the forward and the backward histograms . ^^•^'^ A CuBe 
piston-cylinder pressure cell was used with Dafne oil as 
a pressure transmitting medium. The maximum pres- 
sure achieved was 1.4 GPa at 3 K. The pressure was 
measured by tracking the superconducting transition of 
a very small indium plate used as a manometer (cali- 
bration constant for In: dT^/dP = -0.364 K/GPa). In 
order to avoid charge transfer effects due to chain reorder- 
ing in pressurized YBa2Cu30a;, the samples were cooled 
down below 100 K for the /.tSR measurements within less 
than 1 hour after application of the pressure. This time 
is much shorter than the time constant r = 27.7 h (at 
room temperature) for the pressure activated chain re- 
ordering process4^ Below 100 K r is much longer than 
the typical measurement time of a sample (< 24 h)4^ 

High energy muons (p^ ~ 100 McV/c) were implanted 
in the sample. Forward and backward positron detectors 
with respect to the initial muon polarization were used 
for the measurements of the fiSR asymmetry time spec- 
trum A{t) (see Fig. [8]).— Cylindrically pressed samples 
were loaded into the cylindrical CuBe pressure cell. The 
sample dimensions (diameter 5 mm, height 15 mm) were 
chosen to maximize the filling factor of the pressure cell. 
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FIG. 1: (Color online) /iSR asymmetry signal A{t) of 
YBa2Cu306.98 measured at T = 4.5 K and 95 K in an ap- 
plied field -Bapp = 0.1 T (empty and full circles, respectively). 
The fast relaxation of the fiSK signal (empty circles) is due to 
the formation of a vortex lattice in the superconducting state. 
The solid lines are fits of the data to Eq. fl}. For a better 
visualization the spectra and the fits are shown in a rotating 
reference frame of 0.08 T. 



The fraction of the muons stopping in the sample was 
approximately 40%. 



III. RESULTS AND ANALYSIS DETAILS 



For type-II superconductors in the vortex state in an 
applied field of -Bapp ^ Bc2 {Bc2 is the upper criti- 
cal field) the square root of the second moment of the 
muon depolarization rate a is inversely proportional to 
the square of the magnetic penetration depth: a oc 1/A^ 
(Rcfs. I4l l46ll47l ) and therefore directly related to the su- 
perfluid density: ps oc 1/A^ oc a. For a polycrystalline 
sample of a highly anisotropic and uniaxial superconduc- 
tor the dominant contribution to the muon depolariza- 
tion originates from the in-plane magnetic penetration 
depth Aab = Aeff/1.31, where AeS is an effective (aver- 
aged) magnetic penetration depth4^i^ 

As was pointed above a substantial fraction of the 
/iSR asymmetry signal originates from muons stopping in 
the CuBe material surrounding the sample. The sample 
in the superconducting state induces an inhomogeneous 
field in its vicinity (see Appendix). This leads to an ad- 
ditional depolarization of the /iSR signal arising from the 
muons stopping in the pressure cell. Therefore, the /.iSR 
asymmetry time spectra are characterized by two compo- 



ncnts and may be described by the following expression: 
A{t) =Ai • exp ( -^(^' + '^DtA cos{-f^Bit + 0)+ (1) 



A2- I P{B')cos{jf,B't + (j))dB'. 



Here, Ai and A2 are the initial asymmetries of the two 
components of the /iSR signal {Ai: sample, A2: pres- 
sure cell), 7^ is the gyromagnetic ratio of the muon 
(7^ = 27r X 135.5342 MHz/T), and (j) is the initial phase 
of the muon spin polarization. Bi is the field in the cen- 
ter of the sample (or approximately the mean field in the 
sample). The parameter a denotes the muon depolariza- 
tion in the sample due to the field distribution created by 
the vortex lattice, while (t„ = 0.10(2) /xs~^ is a temper- 
ature, doping, and pressure independent depolarization 
rate due to the nuclear moments present in the sample. 
The total asymmetry is Ai + A2 = 0.275 at 0.1 T and 
0.265 at 0.5 T with Ai/{Ai + A2) ~ 0.4 (~ 40% of the 
muon ensemble are stopping inside the sample). P{B') 
represents the magnetic field distribution probed by the 
muons stopping in the pressure cell as described in detail 
in the Appendix. 

Figure [T] exhibits ^SR asymmetry time spectra of 
YBasCuaOg.gs above (T = 95 K) and below (T = 4.5 K) 
the superconducting transition temperature T^ — 89.6 K 
obtained in an applied field of 0.1 T. For a better visu- 
alization the spectra and the fits are shown in a rotating 
reference frame of 0.08 T. Above Tc only a weak depo- 
larization of the muon spin polarization is visible;^ while 
below Tc the strong relaxation of the /iSR signal reflects 
the formation of the vortex lattice in the superconduct- 
ing statci^i^iii^'^ Figures [2^, b, and c show the Fourier 
transforms (FT) of the /iSR time spectra shown in Fig. 
[TJ In Fig. Hi the FT spectra of YBasCusOe.e below 
and above Tc = 60 K are also shown. The narrow sig- 
nal around Bapp = 0.1 T in Fig. [^b originates from the 
pressure cell, while the broad signal with a first moment 
significantly lower than i?app arises from the supercon- 
ducting sample. It can be seen that the signal of the 
pressure cell is also modified below Tc due to the dia- 
magnctic response of the superconducting sample. The 
solid lines are the FTs of the fits to the data using Eq. ([1]) 
(see also Appendix). The good agreement between the 
fits and the data demonstrates that the model used here 
describes the data rather well. 

The whole temperature dependence of the ^SR asym- 
metry time spectra was fitted globally with the common 
parameters -Bapp, Ai, A2, and an- Solely the parameters 
Bi and a were considered as temperature dependent free 
parameters. As shown in the Appendix the field in the 
sample is macroscopically inhomogeneous due to the in- 
homogeneity of demagnetization effects. Bi is the field 
at the point x = y = z = {i.e., the center of the sam- 
ple). In addition, the parameters describing the muon 
stopping distribution Xo^i and cr, were kept the same for 
each temperature scan (see Eqs. (|A.3[) and (|A.4p in the 
Appendix) . 
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FIG. 2: (Color online) Fourier transform (FT) amplitude as 
a function of field for the spectra shown in Fig. [T] [panels (a) , 
(b), and (c)]. Panel (b) is the expanded [along j/— axis] view of 
panel (a) to show the signal from the sample. Panel (c) is the 
expanded [along a:— axis] view of panel (a) to show the signal 
of the pressure cell. Panel (d) shows the FT of the sample 
with X = 6.6 below and above Tc = 60 K. The solid lines are 
the FTs of the fitted curves shown in Fig. [l] The FT spectra 
are slightly broadened due to a FT apodization of 4 ns~^ . 



The temperature dependence of the depolarization 
rates cr for x = 6.98, 6.8, 6.6, and 6.45 at J5app = 0.1 
and 0.5 T obtained with Eq. ^ are shown in Figs. [3] 
and m respectively. The black empty points correspond 
to the data measured at zero pressure, while the full red 
points correspond to the data measured at 1.1 GPa (for 
X = 6.45, 6.6, and 6.8) and 1.4 GPa (for x = 6.98). 
The values of Tc and cr(0) are in good agreement with 
previous resultsi^'^'^ It is known that the order param- 
eter in YBa2Cu306.98 has predominantly the form of 
A = Ao(p^ ^ Py) [Pi = Pi/\p\ denotes component of the 
unit momentum vector in the reciprocal space along the 
i-th axis] liii^Si^ This implies a linear temperature de- 
pendence of the superfluid density ps down to very low 
temperatures due to quasiparticle excitations at the gap- 
less line nodes in the Px = i \py \ directions on the Fermi 
surface 4^ However, in Fig. [3] we clearly see that <j{T) 
tends to saturate at low temperatures for YBa2Cu306.98 



TABLE I: Summary of the results obtained from the temper- 
ature dependence of a at 0.1 and 0.5 T in YBa2Cu303; using 
Eq. ((2)1. Note that for the sample with x = 6.45 a precise 
analysis of a was not possible due to the occurrence of spin- 
glass magnetism below T ~15 K. Hence, the errors of these 
values of cr(0) are rather large. 
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for both applied magnetic fields. Such a behavior was of- 
ten observed in /iSR studies of polycrystalline samplea^ 
and was explained as originating from a strong scatter- 
ing of electrons on impuritiesi^"— This scattering can 
strongly influence the temperature dependence of ps, but 
it has a minor effect on the superconducting transition 
temperature Tc- In previous theoretical works it was sug- 
gested that such a behavior indicates scattering in the 
unitary limit . ^"^j^^ Thus, the temperature dependence of 
the superfluid density ps was analyzed with the "dirty d- 
wave model" of the BCS theory in the unitary limit of 
carrier scattering as described in Ref. [S^: 
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Here, \ab is the in-plane magnetic penetration depth, 
A(0) = Ao cos(2(^) • g{t) {t = T/T<,) is the 2D-gap- 
function, and e„ — Z(e„)e„ are impurity renormalized 
Matsubara frequencies: e„ ~ {2n + l)'nT. Aq is the max- 
imum of the gap function on the Fermi surface and g{t) 
represents the temperature dependence of the gap with 
5(0) = 1. The parameters Nf and Vf are the density of 
states at the Fermi level and the Fermi velocity, respec- 
tively. The constant e and c represent the electron charge 
and the speed of light. The coefficients Z{en) ares^ 
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FIG. 3: (Color online) (a) Temperature dependence of g of 
YBa2Cu30a: measurcd at Bapp = 0.1 T at zero and applied 
hydrostatic pressures for x = 6.45 (0: -P = GPa; ♦: P = 1.1 
GPa), X = 6.6 (V: P = GPa; T: P = 1.1 GPa), x = 6.8 
(A: P = GPa; k: P = 1.1 GPa), and x = 6.98 (o: P = 
GPa; •: P = 1.4 GPa). The data were analyzed with 
Eq. (lU. The solid curves are fits to the data with Eq. (01. 
(b) Diamagnetic shift of the field /S.B = Bi — -Bapp in the 
corresponding samples. Bi is the mean field in the centrer of 
the sample (see text and Appendix). 
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and Sq ~ 7r/2 in the unitary limit. The angular brackets 
(...) denote averaging over the Fermi surface. In order 

to find Z{en) and g{t), Eq. ([3]) is solved together with 
the following equation>^ 
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Here, i{j{x) is the digamma function. Note that the im- 
purity scattering influences mainly e„ while the temper- 
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FIG. 4: (Color online) Temperature dependence of a of 



YBa2Cu30:E measured at Ba 



0.5 T. The meaning of the 



symbols and the solid lines are the same as in Fig. |3ja 



ature dependence of the gap g(€) changes only slightly 
for a reasonable scattering rate r„. In the clean limit 
(ie., Fti = and Z{en) = 1, Vn) the normalized func- 
tion g{t) is very close to the analytical approximations 
derived from BCS theory.— 

Fits of Eq. © to a(T) oc l/Aab(T)2 measured at var- 
ious hydrostatic pressures are presented in Figs. [3] and 
m The corresponding values for Ag, T^, (Tq, and F^ ob- 
tained from the analysis are summarized in Table HI The 
data for zero and applied pressure and the same doping 
X were analyzed simultaneously with the common pa- 
rameter Tu which characterizes the relaxation rate of the 
Cooper pairs on impurities. As shown in TableUthe data 
for the underdoped samples {x = 6.45, 6.6, and 6.8) are 
well described by the clean limit d-wave model, while for 
the overdoped sample (x = 6.98) F„ = 15(5) K. Here, we 
note that all the studied samples originate from the same 
batch and have an identical thermal history, except of the 
last process of the oxygen reduction. Therefore, we can- 
not explain why only the sample with x = 6.98 exhibits 
a saturation of a in the low temperature limit and why 
it has such a high scattering rate F^ = 15(5) K. Con- 
sequently, we cannot exclude the possibility of a modifi- 
cation of the order parameter in overdoped YBa2Cu302: 
where the pscudogap state gradually vanishes. Such a be- 
havior was also observed previously in optimally doped or 
overdoped polycrystalline samples of YBa2Cu302:i^i^'^ 
However, in single crystal YBa2Cu302; close to optimum 
doping a linear temperature dependence of 1/A^ at low 
temperatures was also reported i^"'^"^ For the sample with 
X = 6.45 only the data above 15 K were analyzed, since 
below 15 K the occurrence of field induced spin-glass 
magnetic order hinders a precise determination of a. 



FIG. 5: (Color online) Tc vs. cr(0) (Uemura plot) at zero and 
applied pressure for YBa2Cu30a; with x — 6.45, 6.6, 6.8, and 
6.98. The solid line is the Uemura line while the dashed line 
is a guide to the eye. The dotted lines represent the pressure 
effect on Tc and (t(0). 
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FIG. 6: (Color online) Relation between Ao and Tc for 
YBa2Cu30a; with x — 6.6, 6.8, and 6.98. The solid line corre- 
sponds to Ao/kBTc — 3 (weak-coupling BSC superconductor: 
Ao/kBTc — 1.76). Both Ao and Ao/kBTc increase with in- 
creasing pressure. 



IV. DISCUSSION 

The main subject of the present study is the pres- 
sure effect on the superconducting gap Aq and the su- 
perfluid density ps oc a. The Uemura relation^, imply- 
ing the linear relation between Tc and ps for underdoped 
cuprate superconductors, was established soon after the 



300 



250 



200 



150 



100 



' 


1 ' . 


o 0GPa,0.1T 
" • 1.1GPa,0.1T 
a 0GPa,0.5T 
■ 1.1GPa,0.5T 


y 


^§^^ 


y^^ 



1.5 



2.0 



2.5 



a(0) (us') 



FIG. 7: (Color online) The gap Ao as a function of a(0) for 
the underdoped samples of YBagCuaOi, with x = 6.6 and 6.8. 
The linear relation between (j(0) and Ao is better fulfilled 
under hydrostatic pressure than the Uemura relation Tc vs 
a(0) and Ao/fcsTc vs. Tc (see Figs. [S]and|51). The line is a 
guide to the eye. 



TABLE II: Values of Op = {dTc/dP)/{d(j /dP) for the un- 
derdoped YBa2Cu302: samples investigated in this work [x = 
6.45, 6.6, and 6.8). 
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25(18) 
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21(6) 



discovery of HTS^ and is one of the important crite- 
ria which a microscopic theory of HTS should explain. 
The Uemura relation for the data summarized in Ta- 
ble U is shown in Fig. [5] As indicated by the dotted 
lines the slope a^ — {dTd dP) / {da / dP) is systemati- 
cally smaller than that suggested by the Uemura line 
with au = dTcjda ~ 40 K//xs~-^. The values of ap 
for the underdoped samples investiganted in this work 
arc summarized in Table [III Note that due to mag- 
netism below ~ 15 K the error of (t(0) for the sam- 
ple with X = 6.45 is rather large. The weighted mean 
value of ftp ~ 23(4) K//is~^ is a factor of ~ 2 smaller 
than ajj ~ 40 (K/xs~^). Such a substantial deviation 
from the Uemura line (with a lower value of ap) was 
also observed by pressure experiments in YBa2Cu408 us- 
ing a magnetization technique.— This is in contrast to 
pressure effect results obtained for the organic supercon- 
ductor K-(BEDT-TTF)2Cu(NCS)2 which follow the Ue- 
mura rclation.-§^ Interestingly, a slope with a factor two 



smaller than that of the Uemura line was also found by 
OIE studies of cuprate superconductors ji^ This suggests 
a strong influence of pressure on the lattice dynamics. 
It is known that the pressure dependence of the super- 
conducting transition temperature is determined by two 
mechanisms: (i) The pressure induced charge transfer to 
Cu02 planes Auh and (ii) the pairing interaction Vcff 
which depends on pressure i^^i^"—'^ 

For the underdoped samples the former mechanism 
dominates (85-90%) the pressure effect on Tci^i^i^ 
Therefore, one can separate the pressure effect on a 
also in two components Act = AcTch + Aay. The first 
term AcTch — {^/ctu){dTc/dP)P follows the Uemura line 
and is mainly due to the charge transfer to the plane. 
The second term Aery ~ (1/ctp - l/au){dTc/dP)P de- 
scribes the increase of the superfluid density solely due 
to a change of the pairing interaction. This increase 
of the superfluid density is equivalent to a decrease of 
the effective mass of the superconducting carriers, since 
Aav/cr = AXy'^/\-^ = -Ato^/to*^ Therefore, the 
pressure-induced change of the effective carrier mass can 
be written as: 

dln{m*y)/dP = -dlii{Xy^)/dP = -(Acry/cr)/AP (6) 

~ [au/ap - l)idTjdP)/T, 

~ a/Tc GPa^^ 

Here, Tc and a are taken at zero pressure and the value 
of (dTc/dP) ~ 4 K/GPa was used. This value is prac- 
tically doping independent in underdoped YBa2Cu30a; 
for 6.45 < a; < 6.8,^ The quantity AAy^ describes the 
change of the superfluid density solely due to a modi- 
fication of the pairing interaction Vgg by pressure. It 
is remarkable to observe the qualitative agreement be- 
tween d\n{\y^)/dP and that found in OIE studies for 
d\nX/d\nMo at different carrier dopings {dluMo is the 
relative change of oxygen mass).— Indeed, Eq. ([6]) pre- 
dicts that the pressure effect on to*\/ strongly increases 
with decreasing Tc- 

Another interesting result is the quite small pressure 
dependence of a in the ovcrdopcd sample with x ~ 6.98, 
which is approximately a factor of ~ 2 weaker than that 
reported from magnetization measurements)^ In Fig. [S] 
the gap magnitudes Aq for the samples with x ~ 6.6, 6.8, 
and 6.98 are plotted as a function of Tc- For the under- 
doped samples {x = 6.6 and 6.8) both Aq and Aq/ZcbTc 
increase upon increasing applied pressure. This suggests 
an increase of the coupling strength with increasing pres- 
sure. This behavior is different from that found for the 
OIE on Aq, where a proportionality between Aq and T^ 
was found, implying a constant ratio of Ao/ksTcr^ In 
the overdoped sample {x = 6.98), Eq. ([2]) suggests a 
small reduction of the coupling strength with increasing 
pressure. However, as was mentioned above, the absence 
of a linear temperature dependence of a at low temper- 
atures for the sample with x = 6.98 might also indicate 
that the superconducting order parameter is not of purely 
d-wave character i^2i^ This, on the other hand, may in- 
fluence the result for Aq and its pressure dependence. 



In Fig. [7]for the undcrdopcd samples {x = 6.6 and 6.8) 
Aq is plotted vs. cr(0), showing a linear correlation be- 
tween the two quantities. Note, that this correlation does 
not change with the application of hydrostatic pressure. 
This is in contrast to what is observed for the Uemura 
relation Tc vs. o'(O) and Ao/fcsTc vs. Tc (see Figs. [5]and 

ED- 



V. CONCLUSIONS 

The pressure dependence of the magnetic penetration 
depth A of polycrystalline YBa2Cu30a; {x = 6.45, 6.6, 
6.8, and 6.98) was studied by /iSR. The pressure de- 
pendence of the superfluid density ps cc a oc 1/A^ as 
a function of the superconducting transition Tc temper- 
ature does not follow the well-known Uemura relationi^ 
The ratio Up = {dTjdP)/{da/dP) ~ 23(4) K//xs-i is a 
factor of ~ 2 smaller than that of the Uemura relation 
observed for underdoped samples. However, the value of 
Up is quite close to that found in OIE studies)^ indicat- 
ing a strong influence of pressure on the lattice degrees 
of freedom. We conclude that the contribution of car- 
rier doping to the pressure dependence of A is similar to 
the OIE on A. A weak pressure dependence of the su- 
perfluid density ps was found in the overdoped sample 
{x = 6.98). The superconducting gap Aq and the BCS 
ratio Ao/fc^T'c both increase with increasing applied hy- 
drostatic pressure in the underdoped samples, implying 
an increase of the coupling strength with pressure. Al- 
though the Uemura relation does not hold and the BCS 
ratio is increasing with pressure in underdoped samples, 
the relation between Aq and the fiSK relaxation rate a 
is invariant under pressure. Finally, a model to analyze 
TF pSK spectra of magnetic/diamagnetic samples loaded 
into a pressure cell was developed and successfully used 
in this paper (see Appendix), resulting in a substantial 
reduction of the systematic errors in the data analysis. 
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Appendix: Field distribution in a pressure cell 
loaded with a sample with a non-zero magnetization 

Samples with a strong magnetization placed in a pres- 
sure cell with an applied magnetic field induce a mag- 
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FIG. 8: (Color online) (a) Schematic sketch of the /^SR pres- 
sure instrument GPD at the Paul Scherrer Institute: Cylin- 
drical sample (blue); pressure cell (yellow); muon stopping 
distribution (red ellipse) , and forward and backward positron 
detectors (black), (b) Illustration of the surface current on 
a slice of a homogeneously magnetized cylindrical sample. 
The magnetic field induced by this slice is equivalent to the 
magnetic field of the surface currents, (c) Cross section of 
the cylindrical sample and the surface current distribution in 
a;z— plane, (d) Magnetic field map of the surface currents as 
illustrated in panels (b) and (c). 



netic field in the space around the sample. Typical ex- 
amples of such samples are superconductors (strong dia- 
magnets), superparamagnets, and ferro- or ferrimagnets. 
Thus, muons stopping in a pressure cell (PC) contain- 
ing the sample will undergo precession in the vector sum 
of the applied field and the field induced by the sample. 
This spatially inhomogeneous field leads to an additional 
depolarization of the muon spin polarization which de- 
pends on the applied field and the induced field together 
with the spatial stopping distribution of the muons. 

Consider the most simplest case of a sample with the 
shape of a round cylinder of hight H and radius R placed 
into a cylindrical pressure cell with the same internal ra- 
dius R (Fig.IHK). Typical pressure cell radii used for ^SR 
studies are R = 2.5 - 4 mm. In standard transverse field 
(TF) pSK experiments the pressure cell is placed with the 
cylinder axis oriented vertically while the magnetic field 
is applied perpendicular to the cylinder axis of the pres- 
sure cell and the muon beam direction (see Fig. [5]). Let 
us introduce a cartesian coordinate system with the y- 
axis along the sample cylinder axis, and the z-axis along 
the direction of the applied field. Thus, the x-axis is 
along the initial muon beam direction which is perpen- 
dicular to the forward and backward detector planes (see 
Fig. [8|). The origin of the coordinate system is located in 
the center of the sample. 

In an applied magnetic field H (along the z-direction) 



the sample has a magnetization M. This magnetization 
is the source of an induced field H' (r ) . Let us assume that 
H' is much weaker that the applied field H which is the 
case for superconductors in a magnetic field oi ^qH 3> 
Bel {Bel is the first critical field). Thus, one can neglect 
the spatial variation of the magnetization due to the ad- 
ditional induced field: M = M{H + H'{r)) ~ M{H). 
Typically half (or even more) of all the muons arc stop- 
ping in the PC outside of the sample volume. The muons 
stopping in the macroscopically inhomogeneous field of 
the PC contribute to an additional relaxation of the /iSR 
signal. In order to describe the total fiSK time spectrum 
(sample and PC) one has to model the field distribu- 
tion H'(r). For an applied field H :^ H'{y) one can ne- 
glect the influence of H'^ir) and Hy{^) on the /iSR time 
spectrum, since only the z-componcnt H'^ (r) contributes 
significantly to the muon depolarization. The induced 
magnetic field H'(r) created by a cylindrical sample can 
be calculated as follows)^ 



H'(r) 



1 

47r 



3(M- (r-r'))(r-r') 



M 



r-'|5 



r — r' 



'13 



dr' 



(A.l) 
Here, the integral is taken over the sample volume V. 
For a sample with a constant magnetization the three- 
dimensional integral can be replaced by surface integrals. 
Let us take one slice of width dz out of the sample cylin- 
der and divide it into many small squares dA = dxdy (see 
Fig. [8}d). The field created by the elementary cell of vol- 
ume dV = dxdydz with magnetization M is equivalent 
to the field created by the current Iz = Mdz circulating 
within this square slice as shown in Fig. [5^. It is obvious 
that integration of this field over the whole slice volume 
will leave only a current Iz flowing over the perimeter of 
the slice. The total field of the cylinder is the integral of 
the fields created by these slices with constant current Iz 
(see Figs. [8)3 and c). 

According to the law of Bio-Savart the field in a point 
r created by the elementary currents Idi at the surface 
of the cylinder (with coordinates Tg) is)^ 



H'(r) = 



/ [d4 X (r - r,)] 
47r Ir — r=P 



(A.2) 



The integration is taken over the surface S of the sam- 
ple and d(.s is the elementary length on the surface with 
its direction along the current (the subscript s denotes 
quantities related to the surfaces of the sample. 

The spacial magnetic field distribution around the 
ferro/paramagnetic sample calculated with Eq. (jA.2[) in 
x-z plane is shown in Fig. [8]i. The total field in the 
pressure cell is the vector sum of this field and the homo- 
geneous external field. It is obvious from the figure that 
the field along the z-axis is higher (lower) than the exter- 
nal field in a ferromagnet(diamagnet). Along the a;-axis, 
on the other hand, the field is lower(highcr) than the ex- 
ternal field in a ferromagnet(diamagnct). The maximal 
(minimal) induced field in the PC are just on the border 
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FIG. 9: (Color online) (a) Contour plot of the field distribu- 
tion H'-{y,z) in the j/z— plane for a cylindrical sample with 
R = 2.5 mm and _ff = 15 mm (the geometry of the sample 
used in the experiment), (b) Contour plot of the muon stop- 
ping distribution in the yz— plane. The gray area on the top 
of the sample corresponds to the empty pressure cell space 
where no muons stop (this space is filled with a low-density 
pressure transmission medium). The dashed line indicates 
the sample space, (c) Magnetic field profile of H'^ along the 
J/— axis and (d) magnetic field profile of H'^ along the z— axis. 



of the sample/pressure cell along z (x) direction. Note 
that demagnetization effects are naturally accounted for 
by using Eq. (jA.2[) . Since the sample is not elliptical 
this leads to field inhomogenieties within the volume of 
the sample (see Fig. [9]). As an example Fig. |9] shows the 
magnetic field distribution in the yz— plane for a cylin- 
drical sample with iJ = 15 mm and radius R = 2.5 
mm, together with fields along z— and ?/— axes calcu- 
lated with Eq. (|A.2p . Due to demagnetization effects the 
magnetic field profiles within the sample has peaks at 
the top and bottom edges of the sample where the de- 
magnetizing fields are minimal (Fig. [9];). On the other 
hand, the field profile within the sample close to the cen- 
ter is quite homogeneous, since a cylinder with infinite 
hight H is equivalent to an ellipsoid in which the field is 



homogeneous. 

In order to calculate the probability field distribution 
of a sample in a PC with a substantial first moment a 
model for the muon stopping distribution is required. 
This distribution may be well approximated by a three- 
dimensional Gaussian!^ 



Psixi,X2,X3) 



A 



^ ^ l — l 



^ 1 



exp 



{Xi 



Xo,; 



2(7? 



(A.3) 

where the subscripts i ~ 1,2,3 correspond to x, y, or 
z, respectively. The quantities ccg.i determine the mean 
value of the muon stopping distribution, ai are corre- 
sponding standard deviations, and A is the normalization 
factor. The quantities xq.i, a;o,2, and xo,3 can be deter- 
mined quite accurately before starting the experiment 
by tuning the momentum of the muon beam and vertical 
positioning of the sample. For a sample with nearly the 
same density as the pressure cell a;o,i — 10,2 — 3:^0,3 — 0. 
Simulations of the stopping distribution with the SRIM 
software^ yield Ci = 0.875 mm for copper (the basic 
component of the CuBe pressure cell) and the minimal 
ratio of cra/cri = 3.36. A maximal ratio of o-ij a\ ~ 4 is 



estimated for the muon beam collimated by a 4 x 10 mm 
collimator (this uncertainty is related with the degree of 
muon beam focusing). The parameter G2 is in fact the 
standard deviation of the function representing the con- 
volution of a Gaussian with cr = 173 over the collimator 
profile function along the y-axis. These parameters de- 
fine the fraction of muons stopping in the PC and the 
sample for a given sample geometry. For a known Pg (r) 
one can calculate the magnetic field probability distribu- 
tion P(&) in the pressure cell by solving the integral: 

P(B)^ J P,{v)5{B~iio[H + H',{v)])dv. (A.4) 



Here, 5{x) is the delta function. The integration is taken 
over the volume of the pressure cell. Note that this is not 
simply the probability field distribution in the pressure 
cell, but it is weighted with the muon stopping probabil- 
ity distribution Ps{x,y,z). Fits of P{B) to the experi- 
mental /iSR data are shown in FigO The function P{B) 
describes the experimentally measured fiSK signal rather 
well. 
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